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FOREWORD 


The  interdiction  problem  treated  in  this  report  originated  with  Frederic  C.  Clodius  (retired) 
of  the  Naval  Weapons  Laboratory  (now  Naval  Surface  Warfare  Center,  Dahlgren  Division 
[NSWCDD])  in  1970.  Lack  of  adequate  computing  speed  and  computer  memory  requirements 
then  limited  the  generality  and  the  extensive  computing  requirements  of  the  simulation  models 
for  interdicting  airfields.  Readdressing  the  problems  with  the  advanced  computing  capabilities 
of  today  allows  for  more  generality  in  proposed  models. 

The  author  is  indebted  to  Dr.  Robert  J.  McDevitt  (W25)  and  Dr.  Demetrios  Serakos  (W25) 
for  helpful  discussions.  David  Bozicevich’s  (CX7)  scholarly  editing  of  the  report  is  appreciated. 

This  document  was  reviewed  by  John  T.  Beck,  Head,  Ballistic  Missile  Defense  Systems 
Analysis  Branch,  and  David  S.  Richardson,  Head,  Warfare  Systems  Development  Division. 
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GLOSSARY 

The  numeral  at  the  end  of  each  item  is  the  page  number  where  the  item  is  first  used. 

(a)  -  Refers  to  the  first  of  two  problems  considered,  1 

(b) -Refers  to  the  second  of  two  problems  considered,  2 

(c) -Refers  to  five  different  subdivisions  of  side  [(0,0)- (LF,0)] ,  4 

(d) -Refers  to  five  different  subdivisions  of  side  [(0,0)-(0,WF)],  4 
AQj-Identihes  the  subset  of  Cy  that  contains  aim  points,  4 
CV  j-Identifies  the  rectangular  cells  that  partition  F,  4 

Ix(a,  b)-Incomplete  beta  function  ratio  with  parameters  a,b,  9 

r(x)-Complete  gamma  function  of  x,  10 

7-=(e-d)/2,  8 

A-Level  of  conhdence,  8 

ctda “Standard  deviation  of  deflection  aiming  error  (y- direction),  1 
(jDB“Standard  deviation  of  deflection  ballistic  dispersion  of  ijth  bomb,  1 
ctra -Standard  deviation  of  range  aiming  error  (x-direction),  1 
crRB“Standard  deviation  of  range  ballistic  dispersion  of  ijth  bomb,  1 
A-Lethal  radius  of  a  bomb,  1 
AF  airfield  or  subrectangle,  1 

EAC-Denotes  an  enlarged  AC  cell  of  polygonal  shape,  6 
EAF-Rectangle  enclosing  AF  with  sides  ELR,  EWR,  1 
ELR-Length  of  rectangle  EAF,  ELR  =  LR  +  2A,  1 
EWR-Width  of  rectangle  EAF,  EWR  =  WR  +  2A,  1 
F-Original  specified  held  or  runway  to  be  interdicted,  1 
f-Total  number  of  failures  to  interdict  F  in  m  sorties,  3 
IC-Identihes  the  five  subdivisions  of  (c),  takes  values  1-5,  4 
IL-Integer  part  of  LF /LR,  4 
IW-Integer  part  of  WF /WR,  4 

JC-Identihes  the  five  subdivisions  of  (d),  takes  values  1-5,  4 

LA-Circular  lethal  area,  with  radius  A,  of  a  bomb,  1 

LA2-Segment  of  circle,  with  circle  bounding  the  LA,  6 

LF-Input  length  of  the  held,  1 

LR-Length  of  a  subrectangle,  AF,  1 

LRl-Length  of  an  AC  cell,  4 

m-Total  number  of  sorties=s+f,  3 

N(0,l)-Normal  probability  distribution  with  0  mean,  variance=l,  2 
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NB-Number  of  bombs  dropped  in  a  sortie,  1 
NK-Number  of  sorties  used  in  (a),  6 

NP-The  number  of  chords  approximating  a  quarter  circle  of  an  EAC,  6 
NP2-Denotes  the  number  of  vertices  in  an  EAC,  4NP  +  5,  6 
p-=s/m,  probability  of  interdicting  F  using  m  sorties,  3 
RL-Fractional  part  of  LF /LR,  4 
RW-Fractional  part  of  WF /WR,  4 

s-Total  number  of  successes  to  interdict  F  in  m  sorties,  3 
WF-Input  width  of  the  held,  1 
WR-Width  of  a  subrectangle,  AF,  1 
WRl-Width  of  an  AC  cell,  4 
XA(k)-x-coordinate  of  kth  aim  point,  2 
XH(k)~x-coordinate  of  kth  hit  point,  2 
YA(k)-y-coordinate  of  kth  aim  point,  2 
YH(k)-y-coordinate  of  kth  hit  point,  2 
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I.  INTRODUCTION 

In  this  report  we  revisit  the  work  by  Dr.  Milton  P.  Jarnagin,  Jr.,  described  in  [1],  entitled: 
A  Runway  Interdiction  Program. 

A  rectangular  field  F  of  dimensions  (LF  x  WF)  with  sides  parallel  to  the  X  and  Y  axes 
is  specified  with  its  lower  left-hand  corner  at  (0,0)  (see  Figure  1).  “Sortie”  and  “run”  are 
the  terms  used  to  note  the  dropping  of  NB  identical  bombs  directed  at  F.  Interdiction  of 
F  means  that  it  is  impossible  to  find  an  enclosed  undamaged  subrectangle,  AF,  with  sides 
parallel  to  the  axes  and  of  prespecified  dimensions  (LR  x  WR)  after  a  sortie  or  sorties.  Thus, 
according  to  [1],  F  is  not  interdicted  if  no  bomb’s  lethal  area  penetrates  the  interior  of  a 
potential  AF,  where  the  lethal  area,  LA,  of  a  bomb  is  the  circle  of  radius  A  with  center 
at  the  impact  point  of  the  bomb  (see  Figure  2).  It  is  assumed  everything  is  destroyed 
inside  LA  and  nothing  is  destroyed  by  that  bomb  outside  LA.  On  this  basis,  any  bomb 
landing  in  a  hypothetical  rectangle,  EF,  enclosing  F  symmetrically  with  sides  of  lengths 
ELF  =  LF  +  2  A,  EWF  =  WF  +  2  A  will  result  in  some  damage  to  F  since  its  LA  will 
intersect  F  except  for  a  region  at  the  corners  of  EF.  Any  bomb  landing  outside  EF  cannot 
contribute  to  the  interdiction  of  F.  In  more  detail,  to  determine  whether  F  is  interdicted 
requires  that  any  AF  must  sustain  damage  from  at  least  one  bomb.  Hence,  any  AF  under 
examination  is  assumed  to  be  symmetrically  enclosed  by  a  rectangle,  EAF,  with  dimensions 
(ELR  x  EWR),  where  ELR  =  LR  +  2  A,  and  EWR  =  WR  +  2  A,  so  that  AF  is  damaged  if 
and  only  if  the  bomb’s  hit  point  is  on  or  in  the  interior  of  EAF.  This  model,  as  given  in  [1], 
introduces  a  bias,  since  the  distance  from  the  corner  of  an  EAF  to  the  closest  corner  point 
of  its  AF  is  greater  than  A  (see  Figure  3,  upper  right  hand  corner).  Thus  the  LA  of  a  bomb 
near  a  corner  of  an  EAF  would  not  intersect  AF,  but  the  model  assumes  otherwise.  This 
can  result  in  optimistic  probabilities  of  interdiction.  The  problem  is  eliminated  in  our  model 
as  described  in  the  next  section. 

In  Figure  2  below,  an  example  of  F  not  interdicted  is  shown,  with  NB  =  4,  since  an 
undamaged  AF  exists.  In  Figure  3,  with  NB  =  4,  an  example  of  F  interdicted  is  shown 
since  no  undamaged  AF  can  exist.  In  the  figures,  noting  some  of  the  variables,  we  have 
LF=6,  WF=2,  LR=2,  WR=1,  A=.25.  Figures  1,2,  and  3  show  an  EAF. 

Two  problems  are  addressed,  both  treated  by  Monte  Carlo  techniques: 

(a)  Determine  the  probability  of  interdicting  F  with  associated  confidence  limits  by  det¬ 
onating  a  given  number  of  identical  bombs,  NB,  with  given  aim  points  and  lethal  circular 
area  with  radius  A.  The  hit  point  of  each  bomb,  with  an  associated  aim  point,  is  subject  to 
the  same  normal  random  aiming  error  with  assigned  standard  deviations  in  range  and  de¬ 
flection,  ctra,  ctda,  and  normal  ballistic  dispersions  with  specified  input  standard  deviations, 
orb,  ctdb- 
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(b)  Determine  the  average  number  of  runs  or  sorties  required  to  interdict  F  with  a  pre¬ 
specified  probability  and  associated  confidence  limits. 

We  proceed  by  briefly  summarizing  the  Jarnagin  model  followed  by  a  description  of  our 
Monte  Carlo  model  (Section  II),  which  differs  significantly.  Both  approaches  use  the  same 
classical  statistical  analysis. 

The  hit  points  for  NB  bombs,  at  the  center  of  the  circles  as  shown  in  Figures  2  and  3, 
have  coordinates  (XH(k),  YH(k)),  k  =  1,...,NB.  They  are  obtained  in  [1]  as  follows:  The 
input  aim  points  are  specified  by  their  coordinates  (XA(k),  YA(k)).  Then  two  normal  ran¬ 
dom  variables,  ipx,ipY,  are  generated  from  N(0, 1)  using  a  subroutine,  RNOR,  from  [2],  so 
that  the  aiming  errors  in  range  and  deflection  for  every  bomb  are  given  by  ipxa ra,  Ay  uda  , 
where  ctRa  and  aDA  are  standard  deviation  inputs,  as  noted  above.  The  ballistic  disper¬ 
sion  errors  are  modeled  for  each  bomb  by  choosing  2  NB  normal  random  variables  0x(k), 
#Y(k)from  N(0,1)  using  RNOR;  so  that  range  and  deflection  ballistic  errors  are  given  by 
#x(k)  ctrb;  0y(k)  ctdb>  where  orb  and  orb  are  standard  deviation  inputs.  The  hit  points  are 
given  by 

XH(k)  =  XA(k)+V'x<TRA  +  ^x(k)(7RB  (1) 

YH(k)  =  YA(k)  +  ^Y<TDA  +  ^Y(k)(7DB,  k  =  1,  . . . ,  NB.  (2) 

2  -  - 


1  .5 


0.5 


O 


Figure  1.  Shows  F,  A,  AF,  EAF. 


Figure  2.  F  not  interdicted  (A  =  .25,  NB  =  4). 


XH(k)  =  XA(k)+V’x<TRA  +  0x(kWB 

YH(k)  =  YA(k)  +  ^Y^DA  +  ^Y(k)uDB,  k  =  1,  . . . ,  NB. 


WF=  =  2 


A  =  .25 


WR  =  1 


lf  =  e 
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Figure  3.  F  interdicted. 

It  is  assumed  the  bomb-depositing  aircraft  fly  parallel  to  the  X  axis.  The  option  of  the 
aircraft  flying  at  an  angle  f3  to  the  positive  X  axis  is  dealt  with  by  a  rotation  of  axes.  Assume 
an  orthogonal  £ 77  set  of  axes  with  the  same  origin  as  the  XY  axes  where  the  £  axis  is  rotated 
counterclockwise  through  an  angle  /3  from  the  positive  X  axis.  Then  the  kth  hit  point  with 
coordinates  (£(k),7y(k))  is  given,  for  the  corresponding  aim  point  (£A (k) ,  r]\ (k) ,  by 


£(k)  =  £A(k)  +  ipt  crRA  +  0{(k)  crRB  (3) 

??(k)  =  77a (k)  +  'ipr,  o-DA  +  0„(k)  crDB,  k  =  1,  . . . ,  NB,  (4) 

where  the  2  NB  +  2  normal  random  numbers  9^( k),  6r](k)  are  chosen  in  the  same  way 

as  those  appearing  in  (1),  (2).  Thus,  in  the  (X,  Y)  coordinates  we  have 

XH(k)  =  £(k)cos(/3) -77(k)sin(/3)  (5) 

YH(k)  =  £(k)  sin(/3)  +  77 (k)  cos(/3),  k  =  1,  . . . ,  NB.  (6) 


This  option  is  taken  from  [1]  but  is  not  included  in  our  software.  It  would  be  a  straightforward 
addition,  if  it  were  needed. 

A  dropping  of  NB  bombs  will  be  called  a  sortie  or  a  run.  Each  sortie  is  assumed  to 
have  the  same  aim  points,  but,  in  general,  with  different  bomb  hit  points.  From  [1],  for 
a  sortie,  new  sequences  are  constructed  that  reorder  the  XH(k)  and  YH(k)  in  increasing 
order  of  magnitude.  Then,  a  set  of  elementary  inequalities  are  evaluated  based  on  these  new 
sequences  to  decide  whether  F  has  been  interdicted.  If  the  answer  is  “yes,”  the  sortie  is  called 
a  success,  s;  otherwise,  it  is  called  a  failure,  f.  If  m  sorties  are  run,  then  the  ratio  p  =  s/m 
is  taken  as  an  estimate  of  the  probability  of  interdicting  F  by  a  single  sortie  (based  on  m 
sorties).  Then,  as  described  in  [1],  confidence  limits  associated  with  p  and  m  are  established. 
We  defer  this  part  of  the  discussion  until  Section  III  since  we  will  use  the  same  statistical 
classical  analysis. 

The  Monte  Carlo  procedure  as  just  described  is  very  efficient  and  fast,  but  it  has  several 
drawbacks.  Namely, 
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1.  As  mentioned  in  the  Introduction,  a  bomb  hit  point  near  the  corner  of  an  EAF  would 
still  allow  credit  to  the  enclosed  AF  in  spite  of  the  fact  that  the  LA  for  that  particular 
bomb  would  not  intersect  the  AF.  This  can  lead  to  optimistic  values  for  p. 

2.  No  a  priori  plans  on  optimum  aiming  of  the  bombs  in  a  sortie  are  discussed. 

3.  The  requirement  that  if  any  part  of  a  bomb’s  LA  intersects  an  AF,  no  matter  how  little, 
that  AF  is  considered  damaged  also  contributes  to  optimistic  estimates  for  p. 

II.  OUR  MODEL 

Consider  F  made  up  of  a  rectangle,  sides  parallel  to  the  X  and  Y  axes  with  its  lower  left- 
hand  corner  positioned  at  (0,0).  In  our  model,  F  is  partitioned  into  rectangular  cells  with 
sides  parallel  to  the  X  and  Y  axes.  Five  different  subdivisions  of  side  (0,0)-(LF,0),  not 
necessarily  equal,  are  used  and  referred  to  as  (c),  and  five  different  subdivisions  of  side 
(0,0)-(0,WF)  are  used,  referred  to  as  (d),  to  generate  the  partitioning.  Thus,  a  total  of  25 
different  partitionings  are  considered.  The  rectangular  cells  are  identified  for  a  particular 
partitioning  as  shown  in  Figure  4,  but  in  general,  are  not  of  equal  size.  We  will  describe  (c) 
understanding  that  (d)  uses  the  same  rules  with  appropriate  changes  of  variables  such  as 
WF  for  LF.  Before  discussing  the  five  subdivision  cases  under  (c),  some  additional  notation 
is  introduced.  Let 

IL  =  integer  part  of  LF/LR  (LF  and  LR  defined  in  the  previous  section). 

IW  =  integer  part  of  WF/WR  (WF  and  WR  defined  in  the  previous  section). 

RL  =  LF  —  IL  *  LR,  fractional  part  of  LF /LR. 

RW  =  WF  —  IW  *  WR,  fractional  part  of  WF/WR. 

Cy  identifies  the  rectangular  cells  that  are  made  up  from  the  partitioning  of  F  with 
i  =  1, ...,  IWC,  j  =  1, ...,  ILC  (See  an  example  in  Figure  4.). 

ACy  denotes  the  subset  of  the  Cy  cells  each  containing  an  aim  point  with 
i  =  2, ...,  IWC  —  1,  j  =  2, ...,  ILC  —  1  (See  an  example  in  Figure  6  where  the  aim  points 
are  indicated  with  small  filled  circles  at  the  center  of  each  AC  cell.). 

The  size  of  the  AC  cells  is  (LRlxWRl).  IC=1,2,3,4,5,  identifies  the  five  subdivisions  of 
(c).  JC=1,2,3,4,5,  identifies  the  five  subdivisions  of  (d). 

We  now  describe  some  of  the  IC,  JC  cases,  with  emphasis  on  the  AC  cells.  More  detail 
is  given  by  assigning  numerical  values  to  the  defining  parameters.  The  idea  behind  the 
generation  of  the  AC  cells  is  to  make  them  as  large  as  possible,  thus  increasing  the  possibility 
of  a  success,  but  small  enough  so  that  no  AF  can  exist  if  all  the  AC  cells  of  a  partitioning 
are  interdicted.  We  assume  that  A  <  WR1  <  LR1. 

Let  LF=6,  WF=2,  A=.25.  Initial  value  for  i  and  j  of  ACy  is  2.  Then 
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IC  =  1,  JC  =  1,  IL  >  1,  IW  >  1,  Figure  6 
LR  =  2,  WR  =  1,  IL  =  3,  IW  =  2,  RL  =  0,  RW  =  0,  LR1  =  LR/2  =  1, 

WR1  =  WR/2  =  .5,  ILC  =  2IL  =  6,  IWC  =  2IW  =  4,  NB  =  4(IL  -  1)(IW  -  1)  =  8. 

IC  =  2,  JC  =  2,  IL  >  1,  IW  >  1,  RL  <  LR/2,  RW  <  WR/2,  Figure  7 
LR  =  2.5,  WR  =  .85,  IL  =  2,  IW  =  2,  RL  =  1,  RW  =  .3,  LR1  =  LR/2  =  1.25, 

WR1  =  WR/2  =  .425,  ILC  =  2IL+1  =  5,  IWC  =  2IW+1  =  5,  NB  =  (2IL-1)(2IW-1)  =  9. 

IC  =  3,  JC  =  3,  IL  >  1,  IW  >  1,  RL  >  LR/2,  RW  >  WR/2,  Figure  8 
LR  =  2.25,  WR  =  .75,  IL  =  2,  IW  =  2,  RL  =  1.5,  RW  =  .50,  LR1  =  LR/2  =  1.125, 

WR1  =  WR/2  =  .375,  ILC  =  2(IL  +  1)  =  6,  IWC  =  2(IW  +  1)  =  6,  NB  =  4IL  *  IW  =  16. 

IC  =  4,  JC  =  4,  IL  =  1,  IW  =  1,  RL  >  LR/2,  RW  >  WR/2,  Figure  9 
LR  =  3.5,  WR  =  1.2,  IL  =  1,  IW  =  1,  RL  =  2.5,  RW  =  .8,  LR1  =  LR/2  =  1.75, 

WR1  =  WR/2  =  .60,  ILC  =  4IL  =  4,  IWC  =  4IW  =  4,  NB  =  2IL  *  2IL  =  4. 

IC  =  5,  JC  =  5,  IL  =  1,  IW  =  1,  RL  <  LR/2,  RW  <  WR/2,  Figure  10 
LR  =  4,  WR  =  1.65,  IL  =  1,  IW  =  1,  RL  =  2,  RW  =  .35,  LR1  =  LR  -  RL  =  2, 

WR1  =  WR  -  RW  =  1.30,  ILC  =  2IL  +  1  =  3,  IWC  =  2IW  +  1  =  3,  NB  =  IL  *  IW  =  1. 

IC  =  2,  JC  =  3,  IL  >  1,  IW  >  1,  RL  <  LR/2,  RW  >>  WR/2,  Figure  11 
LR  =  2.5,  WR  =  .75,  IL  =  2,  IW  =  2,  RL  =  1,  RW  =  .50,  LR1  =  LR/2  =  1.25, 

WR1  =  WR/2  =  .375,  ILC  =  2IL+1  =  5,  IWC  =  2(IW+1)  =  6,  NB  =  (2IL-1)(2IW)  =  12. 

Also,  it  should  be  noted  that  there  exist  chances  that  a  bomb  falls  outside  an  EAC  cell 
and  contributes  to  a  success.  Such  events  are  not  recognized  in  our  model,  thus  tending  to 
a  conservative  estimate  for  the  probability  of  the  interdiction  of  F.  An  example  is  shown  in 
Figure  5  below. 


O  1  2  3  4  5  6 


Figure  4.  Notation  for  a  celled  held. 
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l_R  =  6,  LR  =  2,  WR  =  2,  WF5  =  1  ,  A  =  0.25 


Figure  5.  An  interdicted  field  with  an  empty  AC  cell,  03,4. 


Let  LA2  denote  the  smaller  segment  of  a  circle  centered  at  the  origin  with  radius  A  that 
is  bounded  by  the  chord  (— %/3/2  A,  A/2)  —  (-\/3/2  A,  A/2)  and  the  circular  arc  subtended 
by  the  chord.  We  say  an  AC  cell  is  interdicted  if  at  least  LA2  of  a  bomb  intersects  the 
cell.  This  is  ensured  if  every  AC  cell  is  surrounded  by  a  rectangle  with  rounded  corners, 
namely  quarter  circles,  separated  from  its  AC  cell  by  a  constant  distance  of  A/2.  In  our 
model,  a  quarter  circle  is  approximated  by  NP  equal,  in  length,  chords  subtending  the  arc. 
Denote  the  polygonal  figure  by  EAC  (or  EACy),  so  that,  if  a  bomb  falls  on  or  inside  an 
EAC,  the  associated  AC  is  considered  interdicted.  This  polygonal  approximation,  with 
NP2  =  4NP  +  5  vertices,  leads  to  a  slightly  pessimistic  probability  of  success.  Note  also  the 
overlapping  of  the  EAC’s.  Thus,  one  bomb  can  interdict  more  than  one  EAC.  Figures  6-11 
show  the  EAC’s  for  particular  IC,  JC  with  EAC2,2  depicted  in  bold.  Hence,  if  all  EAC’s  are 
interdicted,  then  F  is  interdicted,  and  our  sortie  is  a  success. 

Since  the  EAC’s  are  polygons,  we  use  a  routine  from  [2],  LOCPT,  which  determines  for 
an  input  polygon  PO,  in  the  plane,  whether  a  point  (xo,  yo)  is  outside  PO.  Thus,  each  EAC, 
for  a  particular  set  IC,  JC,  is  examined  by  LOCPT  to  determine  which  of  the  hit  points,  as 
given  by  (XH(k),  YH(k)),  k  =  1,2,  ...,NB,  (See  (1),  (2)),  interdict  one  or  more  of  the  EAC’s. 

An  added  feature  of  our  software  is  that  we  can  specify  a  sequence  of  sorties,  NK,  in  a 
single  run  keeping  track  of  the  damaged  parts  of  F.  Thus,  if  NK=3  and  if,  after  the  first 
sortie,  F  is  interdicted,  then  the  run  is  declared  a  success;  however,  if  F  is  not  interdicted,  a 
second  sortie  is  made,  and  if  F  is  interdicted  then,  a  success  is  declared.  However,  if  again  F 
is  not  interdicted,  a  final  sortie  is  made  before  a  run  is  declared  a  success  or  failure.  In  [1], 
this  option  is  not  available  or  discussed. 

The  statistical  analyses  that  follow  in  the  next  two  sections  are  taken  directly  from  [1], 
They  are  also  given  in  [3]  and  [4], 
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LF  =  6,  1  FR  =  2,  W  R  =  2,  W  R  =  1  ,  A  =  0.25,  I  R  1  =  1  ,  WR1  =  0.5 


Figure  6.  IC=1,  JC=1.  An  example  of  a  partitioning  of  F. 


LF  =  8,  I  FR  =  2.5,  WR  =  2,  WR  =  0.85,  A  =  0.25,  LR  1  =  1  .25,  WR 1  =  0.425 


Figure  7.  IC=2,  JC=2.  An  example  of  a  partitioning  of  F. 


LF  =  6,  LR  =  2.25,  WR  =  2,  WR  =  0.75,  A  =  0.25,  LR 1  =  1.125,  WR  1  =  0.375 


Figure  8.  IC=3,  JC=3.  An  example  of  a  partitioning  of  F. 
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l_l=  =  6 ,  LR  =  3.5,  WR  =  2,  WR  =  1  .2,  7\  =  0.25,  1R1  =  1  .✓’5,  WR1  =  0.5 


r 
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Figure  9.  IC=4,  JC=4.  An  example  of  a  partitioning  of  F. 


I_R  =  6,  1  R  =  4,  WF=  =  2,  WR  =  1  .65,  A  =  0.25,  IR  1  =  2,  WR1  =1.3 


r 

A 

K 

) 

Figure  10.  IC=5,  JC=5.  An  example  of  a  partitioning  of  F. 


LF  =  6,  I  R  =  2.5,  W  f-  =  2,  WR  =  0.75,  /\  =  0.25,  I  R 1  =  1  .25,  WR  1  =  0.375 


Figure  11.  IC=2,  JC=3.  An  example  of  a  partitioning  of  F. 


III.  STATISTICAL  ANALYSIS  FOR  PROBLEM  (a) 

The  input  for  Problem  (a)  is:  LF,  WF,  LR,  WR,  A,  ctra,  oda,  ctrb,  ctdb,  A,  7.  The  last 
two  will  be  defined  below.  Once  the  IC,  JC  partitioning  of  F  is  established,  NB  is  known. 
Then  NB  aim  points  are  generated,  one  at  the  center  of  each  EAC  cell. 
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A  sortie  or  run  means  that  NB  bombs  are  dropped  with  one  aimed  at  each  aim  point 
as  shown  in  Figures  6-11.  Then,  using  a  normal  random  number  generator  RNOR,  from 
[2],  the  four  sigma  inputs  (standard  deviations)  are  used  to  determine  the  actual  hit  points 
of  each  bomb  as  given  by  (1)  and  (2).  For  a  success,  each  EAC  must  be  interdicted  by  at 
least  one  bomb;  otherwise,  a  failure  to  interdict  has  occurred.  The  decision  as  to  whether 
a  particular  polygonal  EAC  has  been  hit  is  determined  by  the  subroutine  LOCPT  from  [2], 
which  reports  if  the  kth  hit  point  is  inside,  on  the  boundary,  or  outside  any  of  the  EAC’s.  In 
total,  m  runs  or  sorties  are  made  where  s  is  equal  to  the  number  of  successes  and  f  is  equal 
to  the  number  of  failures  so  that  m  =  s  +  f. 

The  quantity  p  =  s/m  is  taken  as  an  estimate  of  the  probability  of  interdicting  F.  By 
computing  confidence  limits  on  p,  we  can  establish  a  level  of  confidence  in  p  by  running 
enough  sorties.  A  value  of  7  is  chosen  to  achieve  a  given  level  of  confidence.  For  example,  if 
7  =  0.05,  then  we  are  aiming  for  a  confidence  level  of  1  —  7  =  0.95%.  The  lower  and  upper 
confidence  limits  d,  e,  are  specified  by  the  input  parameter  A  so  that  0  <  p  —  d  <  A,  and 
0  <  e  —  p  <  A,  or  with  1  —  7  confidence,  the  true  but  unknown  probability  p0  is  in  (d,  e). 

In  order  to  determine  d,  e,  with  q  =  1  —  p,  we  will  need  the  two  relations  that  follow; 

m 


i=s 


m 


p1  qm  1  =  Ip(s,  m  -  s  +  1),  [5,  page  975], 


(7) 


^  r  W  1  =  1  -  Ip(s+  1,  m  -  s)  =  Iq(m  —  s,  s  +  1). 
i=0  7  / 


Note  (8)  follows  from  (7)  since 


E  7  W’-'+E  “  V  q- 

i=o  A  /  i=s+l  A  / 


=  1. 


(8) 

(9) 


The  quantity  Ix(a,  b)  is  known  as  the  incomplete  beta  function  ratio.  It  is  evaluated  by  the 
subroutine  ISUBX  from  [2],  [6].  The  expression  on  the  left  of  (7)  represents  the  probability 
of  s  or  more  successes  in  m  trials.  The  expression  on  the  left  of  (8)  represents  the  probability 
of  s  or  fewer  successes  in  m  trials.  The  algorithm  for  d  and  e  proceeds  as  follows: 

If  p  <  .5,  then  sufficient  sorties  are  run  (summing  successes,  s,  and  failures,  f)  with  a 
clean  held  on  each  one  until  Iq_A(m  —  s,s  +  1)  <7/2,  then  if  p  <  A  set  d  =  0,  and 
compute  e  and  exit,  else  continue  the  sorties  until  Ip_^(s,m  —  s  +  1)  <7/2,  then 
compute  d  and  e  and  exit. 

If  p  >  .5,  then  sufficient  sorties  are  run,  (summing  successes,  s,  and  failures,  f)  with  a 
clean  held  on  each  one  until  Ip__\(s,  m  —  s  +  1)  <  7/2,  then  if  p  >  1  —  A  set  e  =  1,  and 
compute  d  and  exit,  else  continue  the  sorties  until  Iq_A(nr  —  s,s  +  1)  <  7/2,  then 
compute  d  and  e  and  exit,  [3],  [4], 


9 


NSWCDD/TR— 11/189 


In  order  to  compute  d,  let  f(p)  =  Ip(s,  m  — s  +  1)  —  7/2  1 .  Then  d=p  is  obtained  when  f(p)=0. 
The  use  of  the  Newton- Raphson  procedure  to  find  d  requires  the  derivative  of  f  with  respect 
to  p,  f'(p),  i.e., 


f'(p)  = 


T(m  +  1) 


-1  qm_s 


>0, 


(10) 


r(s)r(m  -s  + 1) 

where  T(x)  represents  the  complete  gamma  function  of  x,  [5,  page  255].  Start  with  p  =  d0  = 
(m^r)  2>  then  in  general  dn+i  =  dn  —  f(dn)/f'(dn),  n  =  0,1,....  Convergence  occurs  when 
|(dn+i  —  dn)/dn|  <  eps,  where  eps  is  a  small  positive  number,  say  5e-5.  In  order  to  compute 
e,  let  g(p)  =  Iq(m  —  s,  s  +  1)  —  7/ 2.  Then  e=p  is  obtained  when  g(p)=0.  The  use  of  the 
Newton- Raphson  procedure  to  End  e  requires  the  derivative  of  g  with  respect  to  p,  g'(p), 


i.e., 


g'(p)  = 


T(m  +  1) 


r(s  + 1)  r(m  —  s) 


pq 


<  0. 


(11) 


Start  with  p  =  e0  =  (^)  3,  then,  in  general,  en+i  =  en  -  g(en)/g'(en),  n  =  0, 1, ... .  Con¬ 
vergence  occurs  when  |(en+i  —  en)/en|  <  eps,  where  eps  is  a  small  positive  number,  say  5e-5. 
See  Figure  12. 

The  subroutine  INVBETA  is  used  to  find  p  given  Ip.  It  is  stored  in  INVB.FOR.  The 
functions  in  (10)  and  (11)  are  computed  using  the  subroutine  BRCOMP.  The  latter  routine  is 
taken  from  [2],  The  total  Fortran  software  is  contained  in  four  hies:  JARN1.FOR  INVB.FOR 
POLY2X.FOR,  OINB.FOR.  POLY2X  contains  the  subroutine  ILIW  that  is  used  to  partition 
F  and  set  the  parameters  IC,  JC.  OINB.FOR  contains  supporting  software  taken  from  [2], 
In  addition,  JARN1.FOR  contains  code  that  generates  input  data  for  a  MATLAB  routine 
ILIWPIC.M  to  plot  a  figure  such  as  shown  in  Figures  6-11  based  on  the  particular  value  of 
IC  and  that  of  JC. 


Figure  12.  An  example  of  an  Ip  curve  and  an  Iq  curve. 

xNote  p  here  has  a  different  meaning  than  s/m. 

2Note  f'^do)  =  0.  The  curvature  is  positive  for  p  <  do,  so  the  Newton  iterates  converge. 

3Note  g'^eo)  =  0.  The  curvature  is  positive  for  p  >  eo,  so  the  Newton  iterates  converge. 


10 


NSWCDD/TR— 11/189 


IV.  NUMERICAL  RESULTS  FOR  PROBLEM  (a) 

In  this  section,  some  numerical  results  are  given  for  problem  (a)  in  Table  1.  Results  are 
tabulated  for  six  different  IC,  JC  cases.  The  input  common  to  all  cases  is:  LF=6,  WF=2, 
A=.25,  NP=3,  ctra  =  -75,  oda  =  -50,  orb  =  -25,  odb  =  -20,  with  a  confidence  interval  and 
level  set  at  A  =  .05,  and  7  =  .05,  respectively. 

In  the  first  line  of  output,  under  the  m  columns,  we  show  300(3)  and  220(44).  In  the  first 
case,  the  interpretation  is  that  runs  of  100  sorties  3  times,  accumulating  s’s  and  f’s,  were 
required  before  the  appropriate  inequalities  in  the  boxed  remarks  of  the  previous  section 
were  satisfied.  A  lower  confidence  limit  of  d  =  0.083,  and  an  upper  confidence  limit  of 
e  =  0.159,  were  obtained.  In  the  second  case,  44  sequences  of  5  sorties  each  were  run  before 
the  appropriate  inequalities  of  the  boxed  section  were  satisfied.  A  lower  confidence  limit  of 
d  =  0.075,  and  an  upper  confidence  limit  of  e  =  0.163  were  found.  In  the  second  line  of 
output,  5000(1)  indicates  a  single  set  of  5000  sorties  was  sufficient  to  satisfy  the  appropriate 
inequalities. 
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Table  1.  Various  Cases  for  Problem  (a) 


Confidence  Limits  d,  e,  and  Probability  of  Interdiction,  p 


LR=2,  WR=1,  IC=1,  JC=1,  NB=8,  Figure  6 


NK 

m 

d 

P 

e 

NK 

m 

d 

P 

e 

1 

300(3) 

.083 

.117 

.159 

1 

220(44) 

.075 

.114 

.163 

1 

5000(1) 

.099 

.108 

.117 

1 

10000(1) 

.102 

.108 

.115 

3 

450(3) 

.639 

.684 

.727 

3 

370(74) 

.639 

.689 

.736 

3  5000(1)  .656  .669  .682 


3  10000(1)  .664  .673  .682 


5000(1) 

.106 

.115 

.124 

1 

10000(1) 

.105 

111 

400(4) 

.658 

.705 

.749 

3 

360(72) 

.658 

.708 

5000(1) 

.709 

.721 

.734 

3 

10000(1) 

.714 

.723 

LR=2.25,  W 

Tl=.75,  ic=; 

e  8 

1 


01 


74 


1 

5000(1) 

.000 

.047 

.054 

1 

10000(1) 

.000 

.049 

.053 

3 

400(4) 

.632 

.680 

.725 

3 

375(75) 

.630 

.680 

.727 

3 

5000(1) 

.657 

.670 

.683 

3 

10000(1) 

.661 

.670 

.680 

LR=3.50,  WR=1.20,  IC=4,  JC=4,  NB=4,  Figure  9 

1 

400(4) 

.272 

.317 

.366 

1 

370(74) 

.267 

.314 

.363 

1 

5000(1) 

.283 

.295 

.308 

1 

10000(1) 

.288 

.297 

.306 

3 

250(5 

.824 

.872 

.911 

3 

240(48) 

.822 

.871 

.911 

3 

5000(1) 

.871 

.881 

.889 

3 

10000(1) 

.872 

.878 

.885 

LR=4.00,  WR=1.65,  IC=5,  JC=5,  NB=1,  Figure  10 

1 

400(4) 

.637 

.685 

.730 

1 

375(75) 

.633 

.683 

.730 

1 

5000(1) 

.697 

.710 

.722 

1 

10000(1) 

.705 

.714 

.722 

2 

200(2) 

.873 

.920 

.954 

2 

190(38) 

.867 

.916 

.951 

2 

5000(1) 

.910 

.918 

.925 

2 

10000(1) 

.915 

.920 

.925 

LR=2.50,  WR=.75,  IC=2,  JC=3,  NB=12,  Figure  11 

1 

200(2) 

.054 

.090 

.139 

1 

190(38) 

.053 

.089 

.139 

1 

5000(1) 

.071 

.078 

.086 

1 

10000(1) 

.074 

.079 

.085 

3 

400(4) 

.647 

.695 

.740 

3 

370(74) 

.645 

.695 

.741 

3 

5000(1) 

.695 

.708 

.721 

3 

10000(1) 

.705 

.714 

.723 
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V.  STATISTICAL  ANALYSIS  FOR  PROBLEM  (b) 

In  this  section  we  consider  problem  (b),  first  stated  in  the  Introduction,  namely,  determine 
the  expected  number  of  sorties  to  interdict  F  with  associated  confidence  limits.  The  idea  is 
to  generate  a  set  of  sequences,  each  sequence  made  up  of  a  number  of  sorties  to  interdict  F. 
All  sorties  have  the  same  aim  points.  For  example,  starting  with  a  clean  target  F,  NB  bomb 
hit  points,  determine  if  F  is  interdicted;  if  not,  then  keeping  track  of  the  hit  points,  another 
sortie  is  made  with  the  same  aim  points,  but  with  different  hit  points.  This  is  continued 
until  F  is  interdicted.  Say  eight  sorties  are  required,  then  Ni  =  8 .  A  clean  target  is  used  to 
repeat  the  above  procedure  with  the  ith  sequence  yielding  a  value  of  Nj.  As  a  result,  a  mean 
and  standard  deviation  are  generated  from  k  sequences,  that  is 

k 

Mk  =  ^Ni/k,  (12) 

i=l 

-  Mk)2 

The  value  of  k  is  determined  by  using  the  Student  t  test  to  conclude  when  a  prespecified 
confidence  level  is  achieved,  [3],  [4],  In  [1],  two  parameters,  7  and  A,  are  introduced  for 
this  purpose.  The  number  1  —  7  specifies  the  level  of  confidence  that  the  true  expected 
number  of  sorties  to  interdict,  Mo,  is  in  the  interval  (Mk  —  A,  Mk  +  A).  For  example,  if 
7  =  .05,  A  =  2,  and  Mk  =  4.2,  then  we  have  that  1  —  7,  implying  a  95%  confidence  that 
the  true  mean,  Mo,  is  in  the  interval  (4.2  —  2,4.2  +  2).  In  applying  the  t  test,  the  following 
quantities  are  used: 

t2  =  (kA2)/cr2,  (14) 

x0  =  (k  —  1)/ (k  —  1+  to),  (15) 

Po  =  1  —  Ix[(k  —  l)/2, 1/2].  (16) 

The  quantity  po  represents  the  probability  that  t  of  the  Student  t  distribution  satisfies  the 
inequality  t  <  to  =  A-v/k/ck-  The  derivation  of  these  equations  is  discussed  in  [1],  [7,  see 
page  4],  The  numerical  procedure,  using  (14)-(16),  with  the  zero  subscript  dropped,  is  as 
follows:  the  number  of  sequences,  k,  is  run  a  minimum  number  of  times,  say  ten;  thereafter, 
k  is  increased  by  some  constant  increment,  say  one,  until  p  >  1  —  7 .  At  that  point,  x  is 
found  by  using  Newton-Raphson  to  solve  the  equation 

_ l~Ix[(k  -  l)/2, 1/2]  =  1  -  7.4  (17) 

4Figure  13  shows  the  graphical  relationship  of  (17).  So  that,  starting  with  (1  +  x)/2,  with  x  from  (15),  convergence  of 
Newton-Raphson  is  assured,  since  <//(x)  and  from  (18)  are  shown  in  [1]  to  always  be  negative,  k  >  6.  The  subroutine 

INVBE  is  used  to  find  x  from  (18).  It  is  stored  in  INVB.FOR. 
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or 

0(x)  =  (1  -  Ix[(k  -  l)/2, 1/2])  -  (1  -  7)  =  7  -  Ix[(k  -  l)/2, 1/2],  (18) 

Knowing  x,  (15)  is  solved  for  t2  ,  namely 

t2  =  (k  —  1)(1  —  x)/x  (19) 

Then,  solving  (14)  for  A  gives  A' 

A' =  tcrk/\/k  (20) 

Thus,  we  have  arrived  at  the  result  that,  with  an  assigned  7  say  7  =  .05,  the  true  expected 
value  Mq  lies  in  the  interval  (Mk  —  A',  Mk  +  A')  with  a  confidence  level  of  1  —  7  =  95(%). 


VI.  NUMERICAL  RESULTS  FOR  PROBLEM  (b) 

In  this  section  some  numerical  results  are  given  for  problem  (b).  Results  are  tabulated 
for  six  different  IC,  JC  cases.  The  input  common  to  all  cases  is:  LF=6,  WF=2,  A=0.25, 
ora  =  0.75,  ct^a  =  0.50,  ctrb  =  0.25,  odb  =  0.20,  NP  =  3  with  a  confidence  level  set  at  95% 
and  7  =  0.05,  with  a  confidence  interval  parameter  A  =  1.0.  In  the  first  line  of  output  in 
Table  2,  one  sequence  was  run  made  up  of  100  sorties,  which  was  more  than  enough  sorties 
to  satisfy  p  >  1  —  7  and  achieve  the  95%  confidence  level.  The  first  line  of  output  in  the 
table  refers  to  100  sequence  (k=100),  with  Mk  =  3.04  with  A'  =  0.29.  For  the  same  input, 
three  more  results  are  listed  for  the  same  input  with  k  =  1000,  5000, 10000.  The  results 
in  the  table  show  that  it  is  advantageous  to  run  a  very  large  number  of  sequences  rather 
than  start  at  a  small  number,  say  10,  and  increase  by  an  increment  of  one  until  the  above 
inequality  is  satisfied. 
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Table  2.  Various  Cases  for  Problem  (b) 

Confidence  Level  7  =  0.05,  and  Confidence  Interval,  A  =  1.0 

LR=2,  WR=1,  A=0.25,  IC=1,  JC=1,  NB=8,  Figure  6 

k 

Mk 

Mk  +  A' 

Ok 

k 

Mk  —  A' 

Mk 

Mk  +  A' 

Ok 

100 

2.75 

3.04 

3.33 

1.442 

1000 

3.10 

3.19 

3.29 

1.526 

5000 

3.11 

3.15 

3.19 

1.528 

10000 

3.11 

3.14 

3.17 

1.557 

LR=2.5,  WR=.85,  IC=2,  JC=2,  NB=9,  Figure  7 

100 

2.67 

2.97 

3.27 

1.521 

1000 

2.90 

2.99 

3.08 

1.464 

5000 

2.92 

2.96 

2.99 

1.389 

10000 

2.91 

2.94 

2.97 

1.375 

LR=2.25,  WR=.75,  IC=3,  JC=3,  NB=16,  Figure  8 

100 

2.47 

2.75 

3.03 

1.388 

1000 

2.84 

2.92 

3.01 

1.385 

5000 

2.86 

2.90 

2.94 

1.365 

10000 

2,86 

2.89 

2.91 

1.379 

LR=3.50,  WR=1.20,  IC=4,  JC=4,  NB=4, 

Hgure  9 

100 

1.88 

2.11 

2.34 

1.154 

1000 

2.16 

2.23 

2.31 

5000 

2.17 

2.20 

2.24 

1.157 

10000 

2.19 

2.21 

2.23 

LR=4.00,  WR=1.65,  IC=5,  JC=5,  NB=1,  Figure  10 

100 

1,28 

1.45 

1.62 

0.833 

1000 

1.39 

1.44 

1.49 

0.797 

5000 

1.38 

1.41 

1.43 

0.738 

10000 

1.38 

1.39 

1.41 

0.737 

LR=2.50,  WR=.75,  IC=2,  JC=3,  NB=12,  Figure  11 

100 

2.63 

2.87 

3.11 

1.195 

1000 

3.00 

3.08 

3.17 

1.376 

5000 

3.00 

3.03 

3.07 

1.361 

10000 

2.98 

3.01 

3.04 

1.344 

In  [1]  the  data  for  two  cases  are  discussed.  The  aim  points  are  not  specified.  The  input 
data  of  [1],  besides  that  shown  in  Table  3,  are:  LF  =  6000,  WF  =  200,  NB  =  4,  A  =  20,  k  = 
100,  ctra  =  250,  oda  =  125,  orb  =  40,  odb  =  20.  Our  results  for  comparison  are  given  in 
Table  3  (see  Figure  14).  The  results  of  [1]  are  Mk  =  4.73,  A'  =  .670,  crk  =  3.125  in  one  case. 
In  the  second  case  with  the  same  input,  the  reference  lists  Mk  =  4.01,  X'  =  .405,  crk  =  2.042. 
The  results  in  [1]  indicate  that,  on  the  average,  16-19  bombs  are  required  for  interdiction, 
whereas  our  results  would  indicate  21-22  bombs  would  be  required.  However,  at  least  LA2 
of  a  bomb  must  intersect  an  AF  for  its  interdiction  (see  page  5),  whereas  in  the  case  of  [1], 
if  any  of  the  LA  intersects  AF,  no  matter  how  little,  interdiction  of  that  AF  occurs.  Thus,  a 
more  appropriate  input  A  for  comparison  would  be,  in  our  case,  to  set  A=30  instead  of  20. 
From  Table  3  one  would  conclude,  on  the  average,  around  19  bombs  would  be  required  by 
our  results. 


15 


NSWCDD/TR— 11/189 


Table  3.  Various  Cases  for  Problem  (b) 

(Comparison  of  Results  with  Results  in  [1]) 

Confidence  Level  7  =  0.05,  and  Confidence  Interval,  A  =  1.0 

LR=4200,  WR=50,  A=30,  IC=5,  JC=5,  NB=6,  Figure  14 

k 

Mk  —  A' 

Mk 

Mk  +  A' 

O’k 

k 

Mk  —  A' 

Mk 

Mk  +  A' 

Ok 

100 

2.85 

3.22 

3.59 

1.840 

1000 

3.15 

3.28 

3.40 

1.983 

5000 

3.22 

3.28 

3.33 

1.942 

10000 

3.23 

3.27 

3.30 

1.934 

LR=4200,  WR=50,  A=20,  IC=5,  JC=5,  NB=6,  Figure  10 

100 

3.21 

3.60 

3.99 

1.969 

1000 

3.64 

3.77 

3.90 

2.150 

5000 

3.71 

3.77 

3.83 

2.173 

10000 

3.75 

3.79 

3.83 

2.179 

Finally  a  conjecture  is  offered  that,  in  the  [1]  input  data,  rather  than  RA  =  250,  ctda  = 
125,  orb  =  40,  odb  =  20,  perhaps  the  more  realistic  values,  namely  RA  =  250,  oda  =  40, 
orb  =  125,  odb  =  20  were  intended.  Our  results  in  this  case  are  shown  in  Table  4,  indicating 
a  need  on  the  average  of  only  11-12  bombs  for  interdiction  of  F. 


Table  4.  Various  Cases  for  Problem  (b) 

(Comparison  of  Results  with  Modified  Results  in  [1]) 

Confidence  Level  7  =  0.05,  and  Confidence  Interval,  A  =  1.0 

LR=4200,  WR=50,  A=30,  IC=5,  JC=5,  NB=6,  Figure  14 

k 

Mk  —  A' 

Mk 

Mk  +  A' 

0k 

k 

Mk  —  A' 

Mk 

Mk  +  A' 

°  k 

100 

1.724 

1.88 

2.04 

0.782 

1000 

1.77 

1.83 

1.88 

0.869 

5000 

1.79 

1.81 

1.83 

0.861 

10000 

1.77 

1.78 

1.80 

0.854 

LR=4200,  WR=50,  A=20,  IC=5,  JC=5,  NB=6,  Figure  14 

100 

1.86 

2.01 

2.16 

0.772 

1000 

2.08 

2.14 

2.20 

1.010 

5000 

2.09 

2.11 

2.14 

0.966 

10000 

2.09 

2.11 

2.13 

0.959 
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LF  =  6000,  LR  =  4200,  WF  =  200,  WR  =  50,  A  =  20,  LR1  =  2400,  WR1  =  25 


Figure  14.  IC=5,  C=5.  A  partitioning  of  F  for  data  from  [1],  Table  3. 

VII.  SUMMARY 

The  purpose  of  this  report  is  to  use  Monte  Carlo  procedures  to  produce  software  that 
supplies  statistical  solutions  to  a  military  problem  which  is  beyond  the  desktop  computer 
capabilities  to  utilize  deterministic  analyses. 

A  future  study  could  be  to  adapt  the  developed  software  to  the  present-day  problem 
of  “interdicting”  a  river  with  a  minimum  of  sensors  (replacing  the  bombs  in  the  present 
report).  More  specifically,  determine  the  probability  of  successfully  scanning  the  river  for 
enemy  moving  objects  with  a  minimum  of  appropriately  placed  sensors. 
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